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ABSTRACT 


Variation in the angle a between a pulsar’s rotational and magnetic axes 
would change the torque and spin-down rate. We show that sudden increases 
in a, coincident with glitches, could be responsible for the persistent increases 
in spin-down rate that follow glitches in the Crab pulsar. Moreover, changes 
in a at a rate similar to that inferred for the Crab pulsar account naturally 
for the very low braking index of the Vela pulsar. If a increases with time, all 
pulsar ages obtained from the conventional braking model are underestimates. 
Decoupling of the neutron star liquid interior from the external torque cannot 
account for Vela’s low braking index. Variations in the Crab’s pulse prohle due 
to changes in a might be measurable. 


Subject headings: dense matter — magnetic helds — stars: magnetic helds 
stars: neutron — pulsars: individual (Crab, Vela) 
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1. Introduction 


Although the rotational behavior of pulsars has been monitored in detail for decades, 
some aspects of how an isolated pulsar spins down are unresolved. Questions surround the 
manner in which angular momentum is removed from the system and how the different 
components of the neutron star interior couple to one another. In particular, most pulsars 
do not slow in a regular fashion, but undergo variations in their spin rates in the form of 
glitches and timing noise. Except for such timing irregularities, a pulsar is expected to 
slow down steadily. For example the vacuum dipole model, a theory of rotational energy 
loss to magnetic dipole radiation ( [Goldreich, Pacini, fc Rees 1971|) , predicts oc — 
where is the pulsar rotational velocity. This spin-down law gives a braking index of 
Uobs = = 3. Timing irregularities make meaningful determination of and hence 

Uobs difficult or impossible in most cases. Consequently, the only braking indices available 
until recently were those for three very young, relatively quiet pulsars: the Crab pulsar 
(^obs = 2.51 ± 0.01; P^jyne, Pritchard, fc Smith 1992|) , PSR B1509-58 (uobs = 2.837 ± 0.001; 
[Kaspi, et al. 1994|) , and PSR B0540-69 (uobs = 2.24 ± 0.04; [Boyd et al. 1995|) . These 
braking indices are in significant disagreement with the vacuum dipole model expectation. 
Plasma in the pulsar magnetosphere could carry angular momentum from the star and 
alter the magnetic structure from the dipolar configuration, giving a braking index as low 
as one (see, e.g., [Michel 1973|) . Another possibility is that the magnetic moment of the star 
changes in time, through either a change in the surface field strength (see, e.g., piandfor'd] 
[Applegate, fc Hernquist 1983[ ; [Muslimov fc Page 1995[) or the angle between the magnetic 
and spin axes (see, e.g., [Beskin, Gurevich, fc Istomin 1984 ; [Link, Epstein, fc Baym 1992 ; 
Pandey fc Prasad 1996[). 


In the Crab pulsar, persistent increases in the spin-down rate are observed to accompany 


glitches (see, e.g., [Lyne, Pritchard, fc Smith 1992| and Table 1). This phenomenon can in 
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principle be explained by either a sudden change in the external torque or in the moment 
of inertia acted upon by the torque. In this Letter we suggest that the persistent increases 
in the Crab’s spin-down rate reflect sudden, glitch-induced reorientations of the star’s 
magnetic axis that increase the external torque. 

Recently, Lyne et al. (1996) attempted to separate the glitch activity of the Vela pulsar 
from its underlying spindown, and obtained a surprisingly small value of If this small 
reflects the true spin-down of the star, the braking index is only 1.4 ± 0.2. This braking 
index, far smaller than those for younger pulsars, suggests that the braking index decreases 
with time. This possibility forces reconsideration of traditional determinations of pulsar 
ages that assume constancy of the pulsar magnetic moment. In this Letter we demonstrate 
that shifts in a at a rate similar to that inferred for the Crab account naturally for Vela’s 
low braking index. 


2. The Crab Pulsar 

The Crab pulsar has produced 6 glitches since 1969 (|Lyne, Pritchard, fc Smith 1992|) . 
During each glitch, the spin rate suddenly increases by AD (see column 2 of Table 1). 
Immediately following each glitch, the magnitude of the spin-down rate D is larger than 
before. Most of the excess spin-down rate decays within weeks of the glitch. However the 
long-term behavior, on which we focus here, shows a persistent shift (or offset) AD in 
the spin-down rate that does not decay (see column 3 of Table 1). This type of behavior 
was most striking after the 1989 glitch, in which the persistent shift in spin-down rate 
was especially large (see Fig. 1). Following the glitch, the spin-down rate remained larger 
than that expected from an extrapolation of the preglitch behavior. In all the glitches, the 
shifts in the spin-down rate persisted to the next glitch and were of the same sign. The 
cumulative effect of these offsets amounted to a total increase in the star’s spin-down rate 
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of 0.07% over 23 years. 


Any proposed mechanism for these persistent spin-down rate offsets must conserve 
angular momentum. The evolution of the angular velocity fl{t) of a rigidly rotating pulsar 
under an external torque is 


d{m) 

dt 


= -K 


ext ? 


( 1 ) 


where I is the total moment of inertia. If the star contains a fluid component that is not 
coupled to the external torque, as discussed below, I is less than the star’s total moment 
of inertia. A persistent shift in the spin-down rate could arise through an increase in 
the external torque (|Gullahorn et al. 1977|; Pemiahski fc Proszyhski 1983| ; [Link, Epstein,| 
fc Baym 1992 ) or from a change in I (Paym fc Pines 1971 ; |Alpar fc Pines 1993|) . The 
external torque could vary through either a change in the angle a between the spin and 
magnetic axes, or in the magnitude of the surface magnetic held, B. Gradual growth of 
the surface held has been considered (see, e.g., [Blandford, Applegate, fc Hernquist 1983| ; 
Muslimov fc Page 1995|) , however, no mechanism has been suggested by which the surface 
held could suddenly increase as a consequence of a glitch. Ruderman (1976) has suggested 
that stresses exerted on the inner crust by the superhuid crack the crust and move plates of 
material toward the equator. Because the magnetic held is frozen to the highly-conductive 
crust, such plate tectonics would increase a. Plate tectonic activity could also redistribute 
matter in such a way as to change the orientation of the principle axes of the star. As the 
star relaxes to a new rotational state, a can change, though the direction of change is not 
obvious. 


We now discuss whether the Crab pulsar’s behavior can be understood in terms of 
sudden changes in a. For illustration we use the vacuum dipole model, in which the 
rotational evolution is given by 


G = 


sin^ a 

QcPI 




( 2 ) 
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If the alignment angle a changes by Aa, the spin-down rate of the star shifts by an amonnt 

An 2Aa 


n 


tana 


( 3 ) 


The average rate of change of Afl/fl dne to glitches over 23 years of observations is (|Ly: 


me, 


Pritchard, fc Smith 1992| ) 


To acconnt for these persistent shifts, a must increase at an average rate of 

(d) 1 c: , 1 

' ^ = — = 1.5 X 10"W yr“b 


tan a 


( 4 ) 


( 5 ) 


The growth time is thus 7 x 10^ yr. At this rate of change, the entire angular shift over 
the entire 10^ yr lifetime of the star would be only ~ tage/hi ~ 1/70 rad, considerably less 
than the upper bound of 7r/2. 

Another way for the star’s spin-down rate to increase is through a decrease in the 
effective moment of inertia I on which the external torque acts. Such a change could occur 
through either a structural change of the star, e.g., the star become less oblate, or through 
a decoupling of a portion of the star’s liquid interior from the external torque. The change 
in the spin-down rate, from eq. [ffl would be 


An 

~n 


AI 


( 6 ) 


A decrease in I solely through structural readjustment cannot produce the behavior seen 
in the Crab pulsar. In this case, angular momentum conservation would require the star 
to always spin more rapidly than had the glitch not occurred. The large persistent shifts 
following the 1975 and 1989 glitches, however, eventually caused the star to spin less 
rapidly than had the glitch not occurred (see Fig. 1). In principle, the Crab’s spin-down 
rate offsets could be due to decoupling of a portion of the star’s superfluid interior from the 
external torque (|Alpar fc Pines 1993| ; [Alpar et al. 1996|) , though quantitative agreement of 
this model with the data has yet to be demonstrated. 
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3. The Vela Pulsar 


The high levels of glitch activity and timing noise typically exhibited in older pnlsars 
have prevented the determination of their braking indices. Recently, however, Lyne et al. 
(1996) attempted to separate the glitch activity of the Vela pnlsar from its nnderlying 
spin-down to obtain the average fl. The inferred braking index is 1.4 ± 0.2, signihcantly 
smaller than the braking indices measnred in the yonnger Crab, PSR B1509-58, and 
PSR B0540-69, snggesting evolntion of nobs- Unlike the determinations of the braking 
indices in these yonng pnlsars, which ntilized data intervals that were free of glitches, the 
measnrement of Vela’s riobs employed data containing many glitches. Hence, the valne of 
Vela’s riobs inclndes the effects of any shifts in H that might arise from glitches. Here we 
note that the small braking index for Vela reported by Lyne et al. (1996) conld be dne to 
increases in a of the sort snggested by the Crab’s persistent shifts in spin-down rate. 

Snppose that a pnlsar spins down according to a braking law of the form 


tl = -K{t)9r, 


(7) 


where, for the vacnnm dipole model, K = siv? a / I and n = 3. UK depends on 
time, the observed valne of the braking index, riobs = will differ from n appearing 

in eq. 0. If K is changing at an average rate K, differentiating eq. ||^ in time gives (see, 
e.g., [Michel 199 1|) 

nk k 

n — riobs = —^— = ^tsd — , (8) 

QK R ^ ^ 


where tgd = H/2|H|, the conventional spin-down age, differs from the trne age if K is 
evolving. Hence, if K is increasing with time, nobs will be less than n. In the context of the 
vacnnm dipole model, the growth rate of a reqnired to explain the braking index of Vela is 


(d) _ 1 _ 3 nobs 

tan a Ta 4:tsd 


3.6 X 10 ®rad yr \ 


(9) 








giving a growth time = 3 x 10^ yr. This rate is remarkably close to the average rate of 
orientation shifts inferred from the persistent spin-rate shifts observed in the Crab (see eq. 

The growth time is consistent with Vela’s age, as < tage- 

The changes in a at the average rate given in eq. 0 may occur suddenly, coincident 
with glitches, as we are suggesting for the Crab. Since the average glitch interval in Vela 
is ~ 3 yr, each glitch would then produce a persistent offset in the spin-down rate of 

~ 2 X 10“'^, comparable to the shifts seen in the Crab. These persistent shifts would 
be very difficult to see in Vela; long after a glitch, typically differs from the average 
spin-down rate by ~ 1%, much larger than the expected persistent shifts assuming the 
changes in a occur at glitches. 


In principle. Vela’s lower braking index could be due to a decreasing moment of inertia 
( |Lyne et at. 1996 ). In this case, K/K = —I/I in the vacuum dipole model, and the required 
rate of change is 

^ 3 - Uobs .. 0.8 


/ 

7 


2t 


sd 


0.8 

tsd 


If Uobs remains low over much of the star’s life {^tsd), then the relative change in the 
moment of inertia would be of order unity (see also |lVfichel 1991|) . Explanations based 
on either decoupling of the interior liquid or changes in the star’s structure cannot 
accommodate such drastic decreases in J. In the standard picture of the neutron star 
interior, the only liquid component that can be decoupled is the inner crust neutron 
superfluid, which comprises at most 10% of the star’s total moment of inertia. Furthermore, 
structural changes in I are limited by the deviation from sphericity due to rotation, which 
is A///<10-^ for Vela. 
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4. Discussion 


We have shown that the persistent offsets in the Crab’s spin-down rate following 
glitches conld be dne to sndden glitch-indnced increases in the angle a between the 
rotational and magnetic axes. Moreover, a similar average growth rate of a acconnts for the 
very low braking index of the Vela pnlsar. In the vacnnm dipole model, the characteristic 
growth times of a in these two pulsars are Tq, ~ 7 x 10^ yr for the Crab and Tq, ~ 3 x 10^ yr 
for Vela. 


Lyne and Manchester (1988) found that pulsars older than ~ 10® yr tend to have 
smaller a than younger pulsars. If our suggestion that a grows in young pulsars is correct, 
and the trend found by Lyne and Manchester (1988) is statistically significant, then the 
evolution of a is not monotonic. That is, in young pulsars (tage~10'^ yr) the orientation 
angle a may grow, whereas in older pulsars it becomes smaller. 

It is possible that changes in a are responsible for nobs < 3 for the Crab, PSR B1509-58, 
PSR B0540-69. Since the data intervals used to determine these braking indices contained 
no glitches, to explain these values of nobs the alignment angle a would have to grow between 
glitches. 


The changes in a inferred from the Crab’s persistent offsets might produce measurable 
changes in the pulse profile. For example, a change of a at a glitch could change the 
duration of the line-of-sight’s traverse through the pulse emission cone. The magnitude of 
the associated change in the total pulsed flux is ~ |Aa|/tci/ 2 , where wi /2 is the half-width 
of the emission cone. Taking a = 86°, determined from the radio data ( Rankin 1990|) , Aa 
for the 1989 glitch (which exhibited the largest observed offset) is ~ 3 x 10“® rad (see 
eq. ^). From these numbers we estimate the relative change in pulsed flux to be ~ 1%. 
We estimate a somewhat smaller change, ~ 0.4%, using the value a = 80° obtained by 
Yadigaroglu and Romani (1996) in their analysis of gamma ray and radio data. Such flux 
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changes may be measurable in either x-ray or optical bands. Sekimoto et al. (1995), using 
GINGA, found an upper limit on the change in the pulsed x-ray emission (1-6 keV) across 
the 1989 glitch of < 1.6%. Moreover, Jones, Smith and Nelson (1980), found that the 
pulsed optical emission varied by ~ 1% over 7 yr. 

We thank G. Baym for helpful comments on this manuscript, and A. G. Lyne for 
providing us with data from the 1989 Grab glitch. This work was performed under the 
auspices of the U.S. Department of Energy, and was supported in part by NASA EPSGoR 
Grant 7)^291471, and by IGPP at LANL. 
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Fig. 1.— The behavior of the Crab pulsar surrounding the 1989 glitch at tg ~ 47747 MJD. 
The points show the difference SQ{t) = Q{t) — between the observed spin rate and 

that expected from a model to the preglitch behavior h2o(^)- A shift in the spin-down rate 
of Afl/fl ~ 4 X 10“*^ persisted until the next glitch. By 40d after the glitch, the pulsar was 
spinning slower than had the glitch not occurred. 
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Table 1. Glitches of the Crab Pulsar®^ 


Date 

Glitch size 

10® AD/D 

Persistent spin-down shift 

lO^AD/D 

1969 Sep 

0.40 

0.04 

1975 Feb 

4.4 

2 

1981 

- 

0.1 

1986 Aug 

0.41 

0.2 

1989 Aug 

8.5 

4 

1992 Nov’’ 

- 

- 


®'From Lyne, Prichard, & Smith (1993). 

’’The behavior following the 1992 glitch is not yet published. 







1 0®5a(t)/a„(tg) 







